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ABSTRACT 

We analysed two-dimensional maps of 48 early-type galaxies obtained with the SAURON and 
OASIS integral-field spectrographs using kinemetry, a generalisation of surface photometry 
to the higher order moments of the line-of-sight velocity distribution (LOSVD). The maps 
analysed include: reconstructed image, mean velocity, velocity dispersion, /13 and /14 Gauss- 
?*■ I Hermite moments. Kinemetry is a good method to recognise structures otherwise missed by 

using surface photometry, such as embedded disks and kinematic sub-components. In the 
SAURON sample, we find that 31% of early-type galaxies are single component systems. 91% 
of the multi-components systems have two kinematic subcomponents, the rest having three. 
In addition, 29% of galaxies have kinematically decoupled components, nuclear components 
C^ I with significant kinematic twists. We differentiate between slow and fast rotators using ve- 

locity maps only and find that fast rotating galaxies contain disks with a large range in mass 
fractions to the main body. Specifically, we find that the velocity maps of fast rotators closely 
resemble those of inclined disks, except in the transition regions between kinematic subcom- 
ponents. This deviation is measured with the kinemetric k^/ki ratio, which is large and noisy 
i>j " in slow rotators and about 2% in fast rotators. In terms of E/SO classification, this means that 

I'jjj , 74% of Es and 92% of SOs have components with disk-like kinematics. We suggest that differ- 

j^ ■ ences in k^/ki values for the fast and slow rotators arise from their different intrinsic structure 

which is reflected on the velocity maps. For the majority of fast rotators, the kinematic axial 
ratios are equal to or less than their photometric axial ratios, contrary to what is predicted with 
isotropic Jeans models viewed at different inclinations. The position angles of fast rotators are 
constant, while they vary abruptly in slow rotators. Velocity dispersion maps of face-on galax- 
ies have shapes similar to the distribution of light. Velocity dispersion maps of the edge-on 
fast rotators and all slow rotators show differences which can only be partially explained with 
isotropic models and, in the case of fast rotators, often require additional cold components. We 
constructed local (bin-by-bin) /13 — V/<t and /14 — V/cr diagrams from SAURON observations. 
We confirm the classical anti-correlation of ^13 and V/a, but we also find that /13 is almost zero 
in some objects or even weakly correlated with V/a. The distribution of ^14 for fast and slow 
rotators is mildly positive on average. In general, fast rotators contain flattened components 
characterised by a disk-like rotation. The difference between slow and fast rotators is traceable 
throughout all moments of the LOSVD, with evidence for different intrinsic shapes and orbital 
contents and, hence, likely different evolutionary paths. 

Key words: galaxies: elliptical and lenticular - galaxies: kinematics and dynamics - galaxies: 
structure - galaxies:evolution 
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1 INTRODUCTION 

The classification of galaxies both acknowledges the complexity 
of these celestial objects and attempts to understand their forma- 
tion and evolution. The Hubble Classification of galaxies (Hubbli 
ll936l ; ISandagelll96lMSandage et al._1975. : .Sandage & Bedke..l9» 
recognises the dichotomy between, broadly speaking, disk and el- 
lliptical galaxies, for historical reasons now often called late- and 
early-types. The classification works well on the late-type galaxies 
in particular, dividing the class into a number of subgroups which 
correlate with properties such as bulge-to-disk ratio, morphology of 
spiral arms, gas and dust content, to name a few, but it fai ls to bring 
a physical insight to our understanding of early-types JTremaind 
[l987), where the classification is based on apparent shape and thus 
dependant on viewing angles. 

In an effort to eliminate the unsatisfactory situation 

iKormendv & Bended ( 11996) proposed a revision of the Hubble clas- 
sification. It was based on two discoveries in the 1970s and 1980s, 
both enabled by an improvement in the techn ical capabilities of as- 
tronomical instruments. A series of papers JBertola & Capacciolil 
1 1 9751 ; IlllingworthI 1 1 977r) showed that bright elliptical galaxies do 
not rotate as fast as they sh ould, if they w ere oblate isotropic sys- 
tems supported by rotation teinnevlll978[) . whereas less bright and 
generally smaller systems, includi ng also bulges of spirals, gen- 
erally agree with such pre dictions iJDavies et alJll983l : lKormendvl 
1 1982 : Kormendv & Illingworth 198^. A co mplementary discov- 
erv feender..l988h ; .Bender et al...l98'a.ll989h that the fast rotating 
galaxies are more likely to have disky isophotes, while the slow ro- 
tating galaxies have boxy isophotes, linked again the kinematics and 
shape of galaxies. Kormendy & Bender ( 1996) changed the unifor- 
mity of early-type galaxies to a dichotomy (disky vs boxy, fast vs 
slow, brighter vs less bright) and linked the whole Hubble sequence 
from right to left, from Sc-Sb-Sa to SO-E types, where 'rotation 
decreases in dynamical importance compared to random motions'. 

This important step forward introduced a readily measurable 
parameter related to some physical properties. However, the higher 
order variations in the isophotal shape (diskiness/boxiness) are not 
meas urable at all inclinat ions regardless of the prominence of the 
disks mix & Whitell99(]|) . and finally, they are used to infer the dy- 
namical state of the galaxy. This might be a decent approximation, 
especially if one assumes that all fast rotating galaxies comprise 
of spheroidal sl ow rotating com ponents and disks seen at differ- 
ent inclination l IRixet alj|l999l) . but the edge-on observations of 
spheroidal componen ts in spiral galaxies showed th at bulges are 
rotating fast as well (Kormendv & lUingworthl 1 1 9821) . To compli- 
cate things further the updated classification continues to distin- 
guish SOs from Es, keepin g a viewing angle dependent definition 
of SOs dvan den Berg"ti|[i99Q) . 

The choice of using the fourth order Fourier term in the isopho- 
tometric analysis for classification is natural, because (i) it is much 
easier to take images of galaxies than to measure their kinematics, 
and (ii) until recently it was not realistic to spectroscopically map 
their two-dimensional structur e. The advent of pa noramic integral- 
field units, such as SAURON ( iBacon et al]|200lh . is changing the 
technical possibilities and the field itself; it is now possible to sys- 
tematically map kinematics of nearby galaxies up to their effec- 
tive radii. We have observed 72 nearby E, SO and Sa galaxies as 
part of the SAURON survey dde Zeeuw et alj|2002l . hereafter Pa- 
per II). Focusing here on a subsample of 48 early-type galaxies (E, 
SO), these observations clearly show the previously hinted rich va- 
riety of kinematic substructures such as: kinematically decoupled 
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cores, kinematic twists , counter-rotating structures and central disks 
JEmsellem et al.ll2004 hereafter Paper III). 

Analysing the global propertie s of the SAU RON velocity and 
velocity dispersion maps, Emsellem et alj ( l2007i hereafter Paper 
IX) were able to separate the early-type galaxies into two physi- 
cally distinct classes of slow and fast rotators, according to their 
specific (projected) angular momentum measured within one effec- 
tive radius, Xr. This finding augments the view that led to the re- 
vision of the classification, but the SAURON observations provide 
the crucial quantitive data. Moreover, the results of Paper IX sug- 
gest a way to dramatically improve on the Hubble classification and 
substitute SOs and (misclassified) disky ellipticals with one class of 
fast rotators. 

ICappellari et alj ( 120071 hereafter Paper X) addressed again the 
issue of orbital anisotr opy of early-ty pe galaxies. They constructed 
the (V/a, t) d iagram teirmevlll978r) using an updated formalism 
JBirmevluOOSh . and compared it with the results from general ax- 
isynmietric dynamical models for a subsample of these galaxies 
JCappellari et alJuOO o. hereafter Paper IV). They found that slow 
and fast rotators are clearly separated on the {V/o, e) diagram (un- 
like Es and SOs), such that slow rotators are round, moderately 
isotropic and are likely to be somewhat triaxial, while fast rotators 
appear flattened, span a larger range of anisotropies, but are mostly 
oblate axisymmetric objects. This finding is in a partial agreement 
with previous studies w hich either found ro und early-type galax- 
ies radially anisotropic jvan der Mare]||l99lh . moderately radially 
anisotropic JGerhard et alj|200l [) or only we akly anisotropic with a 
range of anisotropies for flattened systems JGebhardt et al.ll2003n . 
The results of Paper X, however, clearly show that intrinsically flat- 
ter galaxies tend to be more anisotropic in the meridional plane. 
The models also indicate that the fast rotators are often two compo- 
nent systems, having also a flat and rotating, kinematically distinct, 
disk-like component. 

Dynamical models are often time consuming and difficult to 
construct. Ultimately, one would like to be able to classify galaxies 
by their observable properties only. Is it possible to learn about the 
intrinsic shapes of the early-type galaxies from observations only? 
Surface photometry, being but the zeroth moment of the ultimate 
observable quantity for distant galaxies, the line-of-sight velocity 
distribution (LOSVD), cannot give the final answer. It is necessary 
to look at the other moments of the LOSVD: mean velocity, veloc- 
ity dispersion and higher-order moments, c ommonly parameterised 
by G a uss-Hermite co efficients, h-j, and hn, jvan der Marel & Franxl 
I1993I : lGerharq|l993r) , which measure asymmetric and symmetric 
deviation of the LOSVD from a Gaussian, respectively. 

Indeed, in the last dozen years several studies investigated 
higher moments of the LOSVD of early-ty pe galaxies observ- 
ing them along one or multiple slits ( e.g. 'Bender et al. 1 9941: 
van der Ma rel et al.l 1 19941: Iva n der Marell 1 1994: Kronawitter et alj 

Igood^ 



2OO0I: Koprolin & Zeilingeij^ 



Wegner et al.i i2002i : iHau & Forbes! 120061 : ICorsini et al.l l20oi 
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These studies deepened the dichotomy among early-type galaxies 
showing that fast rotating galaxies with disky isophotes also exhibit 
an anti-correlation between h-j, and V/o. This is consistent with 
these galaxies being made of two components: a bulge and a disk. 
The symmetric deviations, on the other hand are usually smaller 
than asymmetric ones, and somewhat positive in general. In addi- 
tion, the observed higher order moments of the LOSVD can be used 
to constrain the possible merger scenarios of early-typ e galaxies 
and their formation in general (e.g.lBalcellsll99 1 ; Bendo & BarnesI 
I2OOOI : iGonzalez-Garcia et alj|20^ : iNaab et alj|200a : iJesseit et alj 



120070 . However, observations along one or two slits are often not 
able to describe the kinematical rich ness of early-type gala xies. 

In this work we use kinemetry jKrainovic et alj|2006l) . a gen- 
eralisation of surface photometry to all moments of the LOSVD, 
to study SAURON maps of 48 early-type galaxies. The purpose of 
this work is to investigate observational clues from resolved two- 
dimensional kinematics for the origin of the differences between 
the slow and fast rotators. 

In Section[2] we briefly remind the reader of the SAURON ob- 
servations and data reduction. Section[3]describes the methods and 
definitions used in this work. The main results are presented in Sec- 
tion|4] In Section[5]we offer an interpretation of the results and we 
summarise the conclusions in Section|6] In Appendix|A]we discuss 
the influence of seeing on the two-dimensional kinematics and in 
Appendix IbI we present the radial profiles of the kinemetric coeffi- 
cients used in this study. 



2 SAMPLE AND DATA 

In this work we used the data from the SAURON sample which was 
designed to be representative of the galaxy populations in the plane 
of ellipticity, e, versus absolute B band magnitude Mb ■ The sample 
and its selection details are presented in Paper II. In this study we 
focus on the 48 galaxies of the SAURON E-l-SO sample. 

SAURON is an integral-field spectrograph with a field-of-view 
(FoV) of about 33" x 41" and 0'.'94 x 0'.'94 square lenses, mounted 
at the William Herschel Telescope. Complementing the SAURON 
large scale FoV, we probed the nuclear regions of a number of 
galaxies with OASIS, then mounted at Canada-Hawaii-France Tele- 
scope, a high spatial resolution integral-field spectrograph, simi - 
laiiy to SAURON based on the TIGER concept l lBacon et al.ll995h . 
The FoV of OASIS is only 10" x 8", but the spatial scale is 
0"27 X 0'.'27, fully sampling the seeing-limited point spread func- 
tion and providing on average a factor of 2 improvement in spa- 
tial resolution over SAURON. The spectral resolution of OASIS is, 
however, about 20% lower than that of SAURON, and only a sub- 
sample of the SAURON galaxies was observed. 

In this paper we are invest igating the stellar kine matics of 
early-type galaxies. Paper III and lMcDermid et alj yOOg, hereafter 
Paper VIII) discuss the extraction of kinematics and construction 
of maps of the mean velocity V, the velocity dispersion o, and the 
Gauss-Hermite moments h^ and /t4 in great detail. All map s used 
in this work are Voronoi binned JCappellari & CopinluOOBO to the 
same signal-to-noise ratio. The SAURON kinematic data used here 
are of the same kinematic extraction as in Paper X with the latest 
improvement on the template mismatch effects in higher moments 
of the LOSVD. The SAURON mean velocity maps are repeated in 
this paper for the sake of clarity, but we encourage the reader to 
have copies of both Paper III and Paper VIII available for reference 
on other moments of the LOSVD. 



3 METHOD AND DEFINITIONS 

Maps of the moments of the LOSVD offer a wealth of information, 
but also suffer from complexity. It is difficult, if not impossible, to 
show error bars for each bin on the map, and the richness of the 
maps can lead to the useful information being lost in detail. As in 
the case of imaging, it is necessary to extract the useful information 
from the maps to profit from their two-dimensional coverage of the 
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objects. In this section we describe the method used to analyse the 
maps and discuss definitions utihsed throughout the paper. 



3.1 Kinemetry 

iKrainovic et alj ( 120060 presented kinemetry, as a quantitative ap- 
proach to analysis of maps of kinematic moments. Kinemetry is a 
generalisation of surface photometry to the higher-order moments 
of the LOSVD. The moments of the LOSVD have odd or even par- 
ity. The surface brightness (zeroth moment) is even, the mean veloc- 
ity (first moment) is odd, the velocity dispersion (second moment) 
is even, etc. Kinemetry is based on the assumption that for the odd 
moments the profile along the ellipse satisfies a simple cosine law, 
while for the even moments the profile is constant (the same as- 
sumption is also used in surface photometry). Kinemetry derives 
such best-fitting ellipses and analyses the profiles of the moments 
extracted along these by means of harmonic decomposition. It fol- 
lows from this that the application of kinemetry on even moments is 
equivalent to surface photometry resulting in the same coefficients 
for parameterisation of the structures (e.g position angle, ellipticity 
and fourth order harmonics). 

Application of kinemetry on odd maps such as velocity map^ 
provides radial profiles of the kinematic position angle PAiti„, axis 
ratio or flattening, qfci„=b/a (where b and a are lengths of minor and 
major axis, respectively), and odd harmonic terms obtained from 
the Fourier expansion (since velocity is an odd map, even terms are, 
in principle, not present, while in practice are very small and usu- 
ally negligible). In the case of stellar velocity maps, the dominant 
kinemetry term is fci = i/af + h\, representing the velocity am- 
plitude, where ai and b\ are the first sine and cosine terms, respec- 
tively. The deviations from the assumed simple cosine law are given 
by the first higher order term that is not fitted, fcs — \/a\ + 65, usu- 
ally normalised with respect to fci . These four parameters form the 
basis of our analysis because they quantify the kinematical prop- 
erties of the observed galaxies: orientation of the map (a projec- 
tion of the angular momentum), opening angle of the iso-velocity 
contours, the amplitude of the rotation and the deviation from the 
assumed azimuthal variation of the velocity map. For the other mo- 
ments of the LOSVD one could derive similar quantities, depending 
on the parity of the moment. As will be discussed below, we focus 
on kinemetry coefficients that describe velocity maps in detail and 
some specific kinemetry coefficients from the maps of the higher 
order moments. A detailed d escription of the metho d, error analy- 
sis and parameters is given in lKrainovic et alj ( 120061) . 



3.2 Radial profiles 

Kinemetric radial profiles can be obtained along ellipses of differ- 
ent axial ratios and position angles. At each radius there is the best 
fitting ellipse, along which a profile of the kinematic moment will 
have a certain shape: it follows a cosine or it is constant, for odd 



^ It is customary in the literature to refer to the maps of mean veloc- 
ity as velocity fields. Sometimes, due to the specific shape of contours of 
constant velocities, velocity maps are refen'ed to as spider diagrams (e.g. 
Ivan der Kruit & Allenlll97a) . Two-dimensional representations of the next 
moment are, however, usually referred to as velocity dispersion maps. In- 
stead of alternating between fields and maps we choose to call all two- 
dimensional representations of the LOSVD maps: velocity map, velocity 
dispersion map, ^3 map.... 



and even moments, respectively. If this is the case, the higher or- 
der Fourier terms are non-existent or at least negligible for such an 
ellipse. 

In the case of even moments, the best fitting ellipses describe 
the underlying iso-contours, like isophotes in the case of surface 
photometry, or contours of constant velocity dispersion, iso-o" con- 
tours. In the case of odd moments, this is somewhat more difficult 
to visualise, but the axial ratio of the best-fitting ellipse is related to 
the opening angle of the iso-velocity contours: the larger the axial 
ratio, the more open is the spider diagram of the velocity map. 

In this study, kinemetry is used for extraction of parameters in 
the following ways: 

i) We apply kinemetry to SAURON reconstructed images of 
galaxies, which are obtained by summing the spectra along the 
spectral direction at each sky position. This is equivalent to low res- 
olution surface photometry on galaxies from the sample. We focus 
on the photometric position angle VAphot and photometric axial ra- 
tio, related to ellipticity as qpft,ot=l-e. In this case, kinemetry is used 
in its even mode, where even harmonics are fitted to the profiles ex- 
tracted along the best fitting ellipses. 

ii) We use kinemetry to derive radial profiles of the four param- 
eters that describe velocity maps: PAfei„, q^in, k\ and fcs. In this 
case, kinemetry is applied to the maps in its odd mode, when only 
odd Fourier harmonics are fit to the profiles extracted along the best 
fitting ellipses, which are, in general, different from the best fitting 
ellipse of (i). In some cases when it is not possible to determine the 
best-fitting ellipse we run kinemetry on circles (see below). 
Hi) Kinemetry is applied to velocity dispersion maps, using the 
even mode as in (i). In this case, however, the parameters of the 
ellipses used to extract profiles were fixed to the best fitting values 
of surface photometry obtained in (i). 

iv) Maps of Gauss-Hermite coefficients /13 and /14 were also pa- 
rameterised using kinemetry in odd and even mode, respectively. In 
both cases, we used the best fitting ellipses from the lowest odd (ve- 
locity map) and even moment (reconstructed image), respectively. 

Before proceeding it is worth explaining in more detail our 
decision not to use kinemetry to fit the ellipses in some cases. Un- 
der (ii) we mentioned that on some velocity maps it was necessary 
to run kinemetry on circles. In general, the mean stellar velocity 
has an odd parity, and its map, in an inclusive triaxial case, will 
be point-antisymmetric. Certain maps, however, do not follow this 
rule, having no detectable net rotation, e.g. NGC4486, or the in- 
ner part of NGC4550. In the latter case, the zero velocity in the 
inner part can be explained by the superp osition of two cou nter- 
rotating s t ellar c omponents as advocated bv lRubin et al.l ( ll992h and 
IrIx et alj ([1993), where the mass of the counter-rotating compo- 
nent is about 50% of the total mass (Paper X). In other cases the 
non rotation could be a result of dominant box orbits which have 
zero angular momentum. The basic assumption of kinemetry for 
odd kinematic moments therefore breaks down resulting in velocity 
maps that appear noisy and one cannot expect reasonable results. 

In practice, this means that the best-fitting ellipse parameters 
for maps without net rotation will not be robustly determined (de- 
generacy in both PAfci„and qkin) while the higher harmonic terms 
will be large and meaningless. Specifically, fcs will have high val- 
ues. We partially alleviate this degeneracy by first running an un- 
constrained kinemetry fit on stellar velocity maps and identifying 
maps where k^/ki > 0.1 and corresponding radii where it occurs. 
At these radii we re-run kinemetry, but using circles for extraction 
of velocity profiles and Fourier analysis. In this way we set the axial 
ratio qfci„=l in order to break the degeneracy. Although the fcs/fci 
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Figure 1. Left: Voronoi binned stellar velocity dispersion maps of 
NGC2549 (top) and NGC4473 (bottom). Over-plotted lines are isophotes. 
In the lower right corner of each map are values that correspond to min- 
imum and maximum colours on the colour bar. Right: Radial profiles of 
62/12 kinemetry terms for NGC2549 (top) and NGC4473 (bottom). 



term cannot be directly compared with the ks/ki term obtained 
from a best-fitting ellipse, in this case, if there is any indication of 
odd parity in the map, we can still determine the local amplitude of 
rotation fci, and give a good estimate for PAkin- 

The other note refers to items (Hi) and (iv). Although, in princi- 
ple, it would be possible to run kinemetry freely on the velocity dis- 
persion maps, or maps of higher Gauss-Hermite moments (e.g. /13 
and /14), the noise in the data is too high to give trustworthy results 
for the whole sample. By setting the shape of the curve to the best- 
fitting ellipses of the corresponding lowest odd or even moment, 
the harmonic terms of kinemetry quantify the differences between 
these even and odd moments of the LOSVD. 

An example of expected differences can be visualised compar- 
ing the isophotes of the surface brightness and the stellar velocity 
dispersion maps of NGC2549 and NGC4473 (Fig.[T). On both im- 
ages isophotes are aligned with the vertical (y) axis of the maps. In 
the case of NGC2549 the contours of constant velocity dispersion 
seem to be perpendicular to the isophotes, at least within the central 
10", while in the rather unusual case of NGC4473, the high values 
of velocity dispersion have the same orientation as the isophotes. 
The physical explanation of these striking differences should be 
looked for in the internal orbital structure. We postpone this dis- 
cussion to Section |4!4l 

The noise and irregular shape of iso-a contours decrease the 
usefulness of fitting for the velocity dispersion contours. Extracting 
harmonic terms along the isophotes, however, can yield a clear sig- 
nal of the different shape of these two moments. An extracted veloc- 
ity dispersion profile in these two cases (for example along the sec- 
ond brightest isophote shown on Fig. [T) will go through two max- 
ima and two minima. The minima and maxima of these two pro- 
files will be out of phase, because along the major axis in NGC2549 
there is a decrease in velocity dispersion while in NGC4473 there is 
an increase. The decomposition of these two profiles will give dif- 



ferent amplitudes to the harmonic terms. Specifically, 62 (cosine) 
term will be the most in fluenced, because th is term is related to the 
error in the axial ratio | jedrzeiewskj||l987 l). and the shapes such 
as in NGC2549 and NGC4473 will give negative (iso-a rounder 
than isophotes) and positive (iso-tj flatter than isophotes) 62, re- 
spectively. An alternative way to visualise the difference between 
these two maps based on the values of 62 is to consider that a neg- 
ative 62, corresponds to a decrease of the a at the major axis of the 
best fitting ellipse compared to the a measured at the minor axis of 
the best fitting ellipse (NGC2549). In contrast, a positive 62 corre- 
sponds to an increased value of the a at the major axis position of 
the best fitting ellipse (NGC4473). By monitoring these harmonic 
terms it is possible to quantify the shape difference between the ob- 
served zeroth and second moments of the LOSVD. 



3.3 Definition of structures on stellar velocity maps 

A few kinemetric profiles are able to describe a wealth of informa- 
tion from the maps. Specifically, we wish to use them to highlight 
the kinematic structures on the maps and to recognise hidden kine- 
matic components. Here we present a set of quantitative criteria for 
describing features on the stellar velocity maps. Some of the criteria 
are dependent on the quality of the data and they should be modified 
if used on m aps obtained with other IFUs. The following rules were 
presented bv lKrainovic et al.l ( 120061) and Paper IX, but here we list 
them for the sake of clarity. 

A single velocity map can contain a number of kinematic com- 
ponents. Often they are easily recognisable by visual inspection. In 
a quantitive way we differentiate between: 

• Single Component (SC) map: having a radially constant or 
slow varying PAfei„and qfei„profiles. 

• Multiple Component (MC) map: characterised with an abrupt 
change in either: Aq^i^ > 0. 1, or APAfci„ > 10° , or a double-hump 
in k\ with a local minimum in between, or a peak in fcs where 
fcs/fci > 0.02. 

MC maps are clearly more complex than SC maps. The above 
values for changes to the kinemetric coefficients are used to deter- 
mine the extent of each subcomponent (components Ci, C2 and C3 
with radii R\2 and _R23 between them). Each subcomponent can be 
described as being of the following type (limiting values apply for 
the SAURON dataset): 

• Disc-like rotation (DR): defined when the higher order har- 
monic k^/ki < 0.02, while the variation of qfci„ and PAfci„ is less 
than 0.1 and 10°, respectively. Note that this name does not imply 
that the object is a disk intrinsically . 

• Low-level velocity (LV): defined when the maximum of k\ is 
lower than 15 kms"^. A special case is Central low-level velocity 
(CLV) when LV occurs in the central kinematical component on the 
map. 

• Kinematic misalignment (KM): defined when the absolute dif- 
ference between the photometric VAphot and kinemetric position 
PAftin angles is larger than 10°. 

• Kinematic twist (KT): defined by a smooth variation of the 
kinematics position angle VAkin with an amplitude of at least 10° 
within the extent of the kinematic component. 

• Kinematically decoupled component (KDC): if there is an 
abrupt change in VAkin, with a difference larger than 20° between 
two adjacent components, or if there is an outer LV component (in 
which case the measurement of VAki„ is uncertain). A special case 
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of KDCs are Counter rotating cores (CRC) where APAfci„ between 
the components is 180° (within the uncertainties). 

Most of the above definitions are new, arising from two- 
dimensional maps which oiTer a more robust detection of structures. 
The defi nition of KDC i s, however, s imilar to the one used in the 
past (e.g. iBendea l 1988a: Statleg |l99lh , where the two-dimensional 
coverage enables a determination of the orientation of the kinematic 
components. It should be noted that classification of a kinematic 
component as a CLV is strongly dependent on the spatial resolution 
of the instrument. As will be seen later, higher spatial resolution 
can change the appearance and therefore the classification of the 
components. 

Similarly, it should be stressed that the limiting values used for 
these definitions are geared towards the SAURON data. The OASIS 
data, due to different instrumental properties and observing set-up, 
will have somewhat different limiting values, mostly arising in the 
higher order Fourier terms. For example, the mean uncertainty of 
ks/ki term for the OASIS sample is 0.033, significantly higher 
compared to the one for the SAURON sample (0.015). In order to 
treat consistently the two data sets, we adopt a somewhat more con- 
servative value of 0.04 as the limiting values for fcg/fci in definition 
of DR component when estimated from the OASIS data. 

While abrupt changes in the orientation, axial ratio, or velocity 
amplitude are intuitively clear as evidences for separate kinematic 
components, the ks/ki as an indicator of components is more com- 
plex to comprehend. Still, the simple models of two kinematic com- 
ponents rotating at a given relative orientation give rise to k^/ki 
term in the kin emetric expansion in the region where these compo- 
nents overlap (Kr ajnovic et al. 2006). Since we measure luminos- 
ity weighted velocities, the position and extent of the raised ks/ki 
region depends on relative luminosity contributions of the compo- 
nents, marking the transition radii between the components and not 
their start or end. Furthermore, it is also, necessary, to distinguish 
between high fcs/fci due to a super-position of kinematic compo- 
nents (a genuine signal) and high fcs/fci originating from noisy 
maps, such as maps with no detectable rotation (e.g. NGC4486) or 
large bin-to-bin variations (e.g. OASIS map of NGC3379). For our 
data, when the signal in fcs is 10% of fci, we consider the noise too 
high and the k^/ki ratio not usable for detecting individual compo- 
nents. 



3.4 Seeing and quantification of kinematic components 

Robust estimates of the number of sub-components in velocity 
maps and their sizes are influenced by three major factors: data 
quality, physical properties and seeing. While the data quality is de- 
scribed by measurement uncertainties, and in that sense it is quan- 
tifiable to some extent, the other two factors are more complex. By 
'physical properties' we assume physical processes that hide kine- 
matic information from our view, such as specific orientation of the 
object, dust obscuration or simply the fact that we are measuring 
luminosity weighted quantities and we might miss kinematic com- 
ponents made up of stars that constitute a low luminosity fraction 
of the total population. 

The influence of seeing is particularly relevant for subcompo- 
nents in the centres of galaxies. In Appendix |A] we tested the de- 
pendence of the kinemetric coefficients on representative seeings, 
for velocity maps viewed at different orientations. This exercise 
showed that: (i) PAkin and k^/ki are not significantly influenced 
by the seeing, (ii) the amplitude and, to a minor extent, the shape of 
fci are somewhat influenced by the seeing, and (Hi) the axial ratio 



qkin can be strongly influenced by the seeing (Fig. lA2b . In addition 
to these conclusions, the test showed that the inclination of an ob- 
ject is also a factor contributing to the change of the intrinsic qkin, 
and to a minor extent, fci profiles, where higher inclinations are par- 
ticularly influenced by the seeing effects. 

In practice, this means that the change in qkin is a less robust 
indicator of kinematic components. We found that more robust in- 
dicators are abrupt changes in fcs/fci and PAkm profiles, double 
humps in fci profiles or decrease of fci amplitude below our detec- 
tion limit for rotation. We used these as estimates of the sizes of 
kinematic components. It should be, however, noted that the size 
of a component is just a luminosity weighted estimate, originating 
from a super-position of luminosities of individual components, and 
the component can intrinsically extend beyond that radius. Only de- 
tailed dynamical models could give a more robust estimate of the 
internal orbital structure. 



3.5 Determination of global and average values 

In addition to radial profiles we present in this paper a number 
of average quantities. Similar luminosity- weighted quantities have 
already been derived in Papers IX and X: global PAkin, global 
PAphot, average e. In this study we use the velocity maps to 
determine the luminosity weighted average (PAkin), {qkin), and 
(fes/fci) for the whole map and for each kinematic component. 
We also measured the luminosity weighted (PAphot), (qphot) from 
the reconstructed images (both global and for each component), 
(62/00) from velocity dispersion maps and average values of hi 
(measured as the ao harmonic term) from /14 maps. In practice, we 
do this following the expression from Paper IX. The mean (G) of a 
quantity G{R) derived from its sampled radial profiles can be ap- 
proximated with: 

j:^^,q(Rk)F(Rk)G(Rk)(RL 



(G) 



R 



2 ^ 

in.k I 



j:ti<iiR>')nR>.)iRLt.k~R. 



(1) 



in.k' 



where q{Rk) and F{Rk) are the axial ratio and the surface bright- 
ness of the best fit ellipse, with semi-major axis Rk- Eq. (TJ is based 
on an expression defined in Ryden et al. ( 1999). The uncertainties 
of these average values are calculated in the standard way as the 
sum of the quadratic differences between the average value and the 
value at each position Rk . 

In Paper X the glob al PAkin wa s derived using the formalism 
from Appendix C in Kiai novic et al.l (2006). This approach differs 
from the one described here in the sense that it is less sensitive to the 
kinematic structures in the central region, such as abrupt changes 
of PA in case of a KDC. That approach is well suited for making 
global comparisons between PAkin or PAphot, such as global kine- 
matic misalignment, when it is required that they are measured on 
large scales to avoid influence of local perturbations in the nuclear 
regions (e.g. seeing, dust, bars). In this study, however, we want to 
compare the radial properties of different moments of the LOSVD 
and for that reason we use the approach of Paper IX to all mea- 
sured quantities. Note that for the purpose of the direct compari- 
son we measured both kinematic and photometric quantities on the 
SAURON data, in contrast with Papers IX and X. 



4 RESULTS 

In this section we present the results of kinemetric analysis of the 
LOSVD moments maps. We look at the presence of kinematic sub- 
structures in velocity maps (Section [4. It , properties of radial pro- 
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files of PAfcin, fci, G (Section [4.2t . comparison between quin and 
c^phot (Section |43t , the shape difference between isophotes and iso- 
o" contours (Section |4!4l ( and properties of fta and /14 Gauss-Hermite 
moments (Section [4.5t with the purpose to investigate the internal 
structure of SAURON galaxies. Kinemetry probes local character- 
istics of galaxies, and we wish to link those with the global proper- 
ties described in Papers IX and X. In this analysis, the most useful 
is the first moment of the LOS VD, the mean velocity, because it is a 
moment rich in structure and with the strongest signal. We present 
the kinemetric profiles of this moment in Appendix |B] Although 
kinemetry is performed on other moments of the LOSVD, we dis- 
cuss the dominant terms only. 



4.1 Substructures on the velocity maps 

Looking at the kinemetric profiles of 48 SAURON galaxies 
(Fig. IB It the following general conclusions can be made: 

• PAt;i„ profiles are in general smooth and often constant (e.g. 
NGC2974) or mildly varying (e.g. NGC474). In some cases there 
are abrupt changes of up to 180° within 1 - 2" (e.g. NGC3608). 

• Profiles of the axial ratio q^^^ are generally smooth and often 
similar to qj,hot profiles (e.g. NGC1023, NGC3384, NGC4570). 

• There is a variety of k\ profiles, most of them rise and flat- 
ten, but some continue to rise, while some drop (e.g. NGC4278, 
NGC4477 and NGC4546). 

• Considering the ks/ki radial dependence, there are three kinds 
of objects: those that have the ratio below 0.02 along most of the 
radius (e.g. NGC2974), those that have the ratio greater than 0.1 
along most of the radius (e.g. NGC4374) and those that have the 
ratio below 0.02 with one (or more) humps above this value (e.g. 
NGC2549). 

• Objects with fcs/fci > 0.1 along a significant part of the profile 
from SAURON data are all classified as slow rotators in Paper IX 
(e.g. NGC4486). 

A further step in understanding the complex velocity maps can 
be made by applying definitions of kinemetric groups (see Sec- 
tion 13. 3t to the radial profiles. They are summarised in Table [T] 
There are 15 galaxies characterised as SC (31%o the rest being 
MC galaxies (69%) of which 10 have KDC (21%) and 4 CLV (8%). 
Higher resolution observations with OASIS, however, show that all 
SAURON CLVs are in fact small KDCs and, moreover, CRCs (Pa- 
per VIII). This means that there are actually 14 (29%) KDCs in 
the SAURON sample. Kinematic profiles of the OASIS data also 
clearly show structures that are partially resolved in the SAURON 
observations, such as KDC (NGC4621, NGC5198, NGC5982) or 
co-rotating components which often have larger amplitudes of ro- 
tation in the OASIS data, corresponding to the nuclear disks visible 
on the HST images. The effects of specific nuclear kinematics, re- 
lated mean ages of the components and possible different formation 
paths were previously discussed in Paper VIII. 

In addition to this grouping of the velocity maps, we can also 
describe the kinematic components. Out of 15 SC galaxies, 8 are 
DR (53%), 2 are KT (13%) and 4 are LV (27%), while the remain- 
ing galaxy (NGC2768) is a rather special case of a solid body rota- 
tor (see below). The majority of MC galaxies have only two kine- 
matic components (30 or 91%), but there are some with three kine- 



^ Although NGC4550 is made of two counter-rotating and co-spatial disks 
(Section [3.2K it is formally characterised as a SC galaxy due to its low ve- 
locity within the SAURON FoV. 




150 
MAX k, [km/s] 



250 



Figure 2. A relation between the maximum amplitude of rotation and 
luminosity-weighted average of k^/ki ratio, which measures departures 
from the cosine law for velocity proiiles extracted along the best fitting el- 
lipse. All points belong to fast rotators. 



matic components (3 or 9%). Many of the components are similar 
in properties. If the inner component (Ci) is DR then the second 
component (C2) is also a DR. If Ci has a more complex kinematic 
structures (KDC, LV, KT), C2 or C3 will in most cases still be a 
DR. Exceptions are found in a few cases when C2 is not rotating 
and can be described as LV. In some cases components may show 
KT, but their k^/ki ratio is less than or equal to 0.02. Counting all 
galaxies that have at least one component with fcs/fci < 0.02, the 
number of systems with DR-like characteristics rises to 35 (73%). 
In terms of E/SO classification, this means that 74% of Es and 92% 
of SOs have components with disk-like kinematics 

Using the resolving power of the OASIS data we can go even 
further: allowing for large error bars and allowing for considerable 
uncertainty of component boundaries, virtually every galaxy that 
shows rotation (including parts of KDCs) has at least one region 
with ks/ki < 0.02 (0.04 in OASIS). This can be seen on Fig.|2] 
that shows luminosity weighted average of k^/ki ratio versus the 
maximum rotational amplitude fci for all fast rotators measured on 
the SAURON data. The large uncertainties in some cases reflect 
the multiple component nature of fast rotators, bec ause fcs/fci ratio 
rises in the transition region between components l IKrainovic et alj 
I2OO6I) . Notably, both average values and uncertainties rise as the 
maximum rotation velocity decreases. This suggest a more com- 
plex structure (more components, larger difference between com- 
ponents) in galaxies with lower amplitude of rotation. 

We estimated local kinematic misalignment for each kinematic 
component. The results show that a total of 25 galaxies (52% of the 
sample) have some evidence of KM. It should, however, be kept in 
mind that it is difficult to determine the sense of rotation for LV 
components. Ignoring galaxies with a LV in the single component 
or in the second component, we are left with 15 galaxies with kine- 
matic misalignments. Of these, 13 galaxies show misalignment in 
the first component, and 7 in the second component. Moreover, only 
5 galaxies show misalignment in two components, as a global prop- 
erty within the SAURON field of view. 

As it can be seen from the kinemetric profiles there are spe- 
cial cases for which kinemetry is not able to determine the charac- 
teristic parameters robustly. This, in general, occurs when velocity 
drops below ~ 15 kms"^ (e.g. NGC4486, NGC4550), and it is 
not surprising since these maps also often lack odd parity, an ex- 
pected property of the first moment of the LOSVD. As discussed 
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Table 1. Kinemetric properties of the 48 E and SO SAURON galaxies. 



Galaxy Group Nc Ri2 R23 Ci C2 C3 KM 

(1) (2) (3) (4) (5) (6) (7) (8) (9) 



Rotator Comment 
(10) (11) 



0474 


MC 




2 


7 


- DR 


DR 


KM(1,2) 


F 


0524 


SC 




1 


- 


- DR 


- 


- 


F 


0821 


SC 




1 


- 


- DR 


- 


- 


F 


1023 


SC 




1 


- 


- KT 


- 


- 


F 


2549 


MC 




2 


13 


- DR 


DR 


- 


F 


2685 


SC 




1 


- 


- DR 


- 


- 


F 


2695 


MC 




2 


7 


- DR 


DR 


- 


F 


2699 


MC 




2 


6 


- DR 


DR 


KM(1,-) 


F 


2768 


SC 




1 


- 


- 


- 


- 


F 


2974 


SC 




1 


- 


- DR 


- 


- 


F 


3032 


MC (CLV) 


2 


2.5 


- LV 


DR 


- 


F 


3156 


SC 




1 


- 


- DR 


- 


- 


F 


3377 


SC 




1 


- 


- KT 


- 


- 


F 


3379 


SC 




1 


- 


- DR 


- 


- 


F 


3384 


MC 




2 


10 


- DR 


DR 


KM(-,2) 


F 


3414 


MC 


(KDC) 


2 


10 


- DR 


LV 


KM(-,2*) 


S 


3489 


MC 




2 


6 


- DR 


DR 


- 


F 


3608 


MC 


(KDC) 


2 


10 


- DR 


LV 


KM(-,2*) 


S 


4150 


MC (CLV) 


3 


3.5 9.5 LV 


DR KM(1,-,-) 


F 


4262 


MC 




2 


9 


- DR 


- 


KM(1,2) 


F 


4270 


MC 




2 


6 


- 


- 


- 


F 


4278 


MC 




2 


16 


- DR 


- 


KM(-,2) 


F 


4374 


SC 




1 


- 


- LV 


- 


KM* 


S 


4382 


MC 


(CLV) 


3 


2 14.5 LV 


DR DR KM(1,-,-) 


F 


4387 


MC 




2 


7 


- DR 


DR 


- 


F 


4458 


MC 


(KDC) 


2 


3 


- 


LV 


KM(-,2*) 


S 


4459 


MC 




2 


12 


- DR 


DR 


KM(1,-) 


F 


4473 


MC 




2 


10 


- DR 


- 


- 


F 


4477 


SC 




1 


- 


- DR 


- 


KM 


F 


4486 


SC 




1 


- 


- LV 


- 


KM* 


S 


4526 


MC 




2 


11 


- DR 


- 


- 


F 


4546 


MC 




2 


9 


- DR 


DR 


KM(1,-) 


F 


4550 


SC 




1 


- 


- LV 


- 


KM* 


S 


4552 


MC (KDC) 


2 


4 


- KT 


- 


KM(1,2) 


S 


4564 


SC 




1 


- 


- DR 


- 


- 


F 


4570 


MC 




2 


8 


- DR 


DR 


- 


F 


4621 


MC 


(KDC) 


2 


4 


- DR 


DR 


KM(1,-) 


F 


4660 


MC 




2 


7 


- DR 


DR 


- 


F 


5198 


MC 


(KDC) 


2 


2.5 


- 


KT 


KM(1*,2*) 


S 


5308 


MC 




2 


7 


- DR 


DR 


- 


F 


5813 


MC 


(KDC) 


2 


12 


- DR 


LV 


KM(-,2*) 


S 


5831 


MC 


(KDC) 


2 


8 


- DR 


LV 


KM(-,2*) 


S 


5838 


MC 




2 


6 


- DR 


DR 


- 


F 


5845 


MC 




2 


4.5 


- DR 


DR 


- 


F 


5846 


SC 




1 


- 


- LV 


- 


KM* 


S 


5982 


MC (KDC) 


2 


3.5 


- 


LV 


KM(1,-) 


S 


7332 


MC 


(KDC) 


3 


3 


12 KT 


DR DR KM(1,-,-) 


F 


7457 


MC 


(CLV) 


2 


3 


- KT 


DR 


KM(1,2) 


F 



KT between Ci and C2 
Possible C2 beyond r=12'' 
Flat fci profile 
fcs/fci < 0.02, 
C2: Flat fci profile 



C2: Flat fci profile 

r < 10" rigid body rotation; fes/fci < 0.02 for r > 10" 

CRC in OASIS map 

ks/ki < 0.02 over the map 

CRC 

CRC 

KDC in OASIS map with r^ize = 1''5 

Possible there components 
C2: Decreasing fci profile 

CRC in OASIS map, KT between Ci and C2 
Decreasing fci beyond r = 13" 

Ci: possible KT. C'2'- decreasing fci profile. 



C2 : decreasing fci profile 

C2: Flat fci profile 

Two co-spatial counter-rotating disks not detected 

Flat fci = 15 km s~^ over the map 



CRC in OASIS with r~ 1'.'5 

C2:LV between 2.5-10" 
No signature of C2 in fcs/fci 



No signature of C2 in fcs/fci 



C2:LV is between 2-8" 

C2: continuously increasing fci; KDC only in PA change 

C2: continuously increasing fci 



Notes: (1) Galaxy identifier; (2) Kinematic galaxy group: see text for details; (3) Number of kinematic components; (4) Transition radius between the first and 
second components (arcsec); (5) Transition radius between the second and third components (arcsec); (6), (7) and (8) Kinematic group for the first second and 
third components; (9) Local kinematic misalignment between luminosity-weighted averages of PAt,,j„ and PApjjoj: numbers refer to the kinematic component 
and * notes that the PA^^^ was determined in the region with fci mostly below 15kms^^; (10) Rotator class: S -slow rotator, F - fast rotator (11) Comment. 



in iKrainovic et aLl (20061), another exception is the case of solid 
body rotation for which the iso-contours are parallel with the zero- 
velocity curve. In the SAURON sample this is seen in NGC2768. 
Since the velocity iso-contours are parallel the axial-ratio and the 
position angle are poorly defined and in practice strongly influenced 
by noise in the data. Determination of the best-fit ellipse parameters 
for the solid body rotation is degenerate. From these reasons this 



galaxy should also be considered with care when comparing with 
kinemetry results for other objects. 



4.2 Fast vs Slow Rotators 

As stated above, all galaxies with k^/ki > 0.1 along a significant 
part of the profile from SAURON data are classified as slow rota- 
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tors in Paper IX. To some extent this is expected since slow rotators 
in general show very little rotation and kinemetry assumptions are 
not satisfied. The origin of the noise in the maps of slow rotators, 
which generates large fcs/fci, is likely reflecting a special internal 
structure, and, in principal, does not come from technical aspects of 
the observations. Even velocity maps with low amplitude of rota- 
tion could show regular spider diagrams (e.g. disks seen at very low 
inclinations) observed at the same signal-to-noise. The one-to-one 
relation between slow rotators and objects with large higher-order 
harmonic terms is significant since the slow/fast rotator classifica- 
tion is based on both velocity and velocity dispersion maps and re- 
flects the internal structure. 

The OASIS data cover only a small fraction of the effective ra- 
dius and do not show this relationship. The slow rotators NGC3414, 
NGC3608, NGC5813 and NGC5982 show considerable rotation 
because the KDC is covering the full OASIS FoV, while the cen- 
tral regions of the fast rotators NGC2768, NGC3032 and NGC3379 
still have small amplitudes of rotation and high (and noisy) fcs/fci 
ratios. 

The velocity maps of the 12 slow rotating galaxies in our sam- 
ple can be described either as LV or as KDC+LV. In that respect, 4 
slow rotators are SC systems (NGC4374, NGC4486, NGC4550 and 
NGC5846) which do not show any detectable rotation (at SAURON 
resolution) and the other 8 are MC systems where Ci is a KDC and 
C2 is a LV. In between these two cases is NGC4552 with a rather 
constant rotation velocity of 15 kms~^, the boundary level for LV, 
and the Ci between a KDC and a large KT. The other three quarters 
of galaxies in the SAURON sample are fast rotators. Only a third of 
them are described as SC, but as shown above (Section|4!T]l all fast 
rotators have components with kinematics that can be described as 
DR. Moreover, in the case of some slow rotators with small, but not 
negligible rotation in the centres (KDC), the higher resolution OA- 
SIS data were able to ascertain that these components have near to 
DR properties (within often large uncertainties) 

Figure |3] has three panels highlighting most obvious kine- 
matic properties of slow and fast rotators. The top panel shows 
radial variations of the kinematic position angles, PAfei„, which 
are present in various forms, ranging from minor twists in the nu- 
clei, through abrupt jumps at the end of KDCs, to almost random 
changes with radius. However, only PA^in of slow rotators are char- 
acterised by strong and rapid changes. Fast rotators show remark- 
ably constant PAfei„. If twists are present in PAfei„of fast rotators, 
they are small in amplitude (< 30°) and confined to the nuclear 
region in shapes of physically small KDCs (NGC4150, NGC4382, 
NGC4621, NGC7332, NGC7457). Slow rotators on the other hand 
show a much greater amplitude in change of PAfei„. 

It should be stressed again that determination of PA^in for 
slow rotators is much more ill defined than for fast rotators in 
the sense that if there is no rotation, there is also no orientation 
of rotation. The abrupt changes in PAkin are the consequence of 
this in some cases (NGC4374, NGC4486, NGC4550, NGC5846), 
and while one could debate the robustness of measured PAkin, one 
should acknowledge the different nature of these systems from ob- 
jects with a constant PAfei„. 

The difference between the slow and fast rotators is most vis- 
ible in the amplitude of rotation. The middle panel of Figure [3] 
shows the radial profiles of fci kinemetric terms for 48 early-type 
galaxies. Most of the profiles cover up to 1 Re in radius. Slow ro- 
tating galaxies can show a non-zero amplitude of rotation in the 
centres (KDCs), but the amplitude is, in general, not very high and 
towards the edge of the map it is mostly negligible. The only ex- 
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Figure 3. Radial profiles of (from top to bottom) PAki„- (PAfei„}, fci and 
(TO from the SAURON data. (R^fein) are luminosity weighted averages of 
PAj;i„ profiles, ctq profiles are normalised at Re/5. The profiles of slow and 
fast rotators are coloured in red and black, respectively. In the middle panel 
dashed blue lines are over-plotted to guide the eye for the cases with specific 
profiles as mentioned in text. 



10 Davor Krajnovic et al. 



ception is NGC5982, which approaches two fast rotators with the 
slowest rotation in the outer regions (NGC4278 and NGC4473). 

Another characteristic of this plot is the variety of profiles. 
They include: monotonically rising profiles (e.g. NGC3032), pro- 
files with an initial slow rise which turns to a more rapid one 
(e.g. NGC524), a rapid rise to a maximum followed by a plateau 
(e.g. NGC4546), rise to a maximum followed by a decrease (e.g. 
NGC4526), double hump profiles (e.g. NGC4660), flat profiles (e.g. 
NGC821) and, in slow rotators, profiles showing a decrease below 
our detection limit. Keeping in mind that the SAURON sample is 
not a complete sample, among fast rotators there are 17 (47%) with 
increasing profile at the edge of the SAURON map, 9 (25%) have 
flat profiles, 5 (14%) decreasing profiles and 4 galaxies have in- 
termediate (difficult to classify) profiles. Among slow rotators there 
are 3 (25%) galaxies with increasing profiles at the edge of the map, 
the rest being flat and below the detection limit. 

These statistics are influenced by the size of velocity maps 
and the coverage of kinematics components by kinemetric ellipses. 
Clearly, larger scale observations would detect the end of rise in 
amplitude in galaxies that are now observed to have increasing fci. 
Similarly, it is possible that a decrease in k\ could be followed by 
an additional increase or a flat profile at large radii. Still, there are 
two general conclusions for fast rotators: they mostly show increas- 
ing velocity profiles at 1-Re, where the range of maximum velocity 
amplitude spans 200 km s~ ^ . On the other hand, slow rotators have 
velocity amplitude mostly less than 20 kms"^ at IR^. 

The bottom panel of Fig. [3] shows radial velocity dispersion 
profiles, ctq, extracted along the isophotes from the velocity dis- 
persion maps as ao harmonic terms. All profiles are normalised to 
their value at Re/5. This highlights the similar general shape of the 
(To radial profiles. The only outlier is NGC4550 with an ctq profile 
which increases with radius. Most of the other profiles, while dif- 
ferent in detail, show a general trend of increasing ctq towards the 
centre and also have a similar shape. A few profiles are consistent 
with being flat (Aao/AR < 30 kms"^) over the whole profile 
(visible only in fast rotators such as NGC7457). 

If there are any real differences between ao profiles, they are 
apparent for radii smaller than Re/5. There are a few exceptions 
to the general trend: i) profiles with a decrease of more than 5% in 
the normalised ctq within Re/5 (e.g. NGC4382), ii) profiles that are 
flat to within 5% inside the Re/5 (e.g. NGC7457), and Hi) profiles 
with a central rise followed by a drop and consecutive rise forming 
a profile with double maxima (e.g. NGC5813). These cases occur 
mostly in fast rotators, with a few exceptions in slow rotators. 

These central plateaus and drops are interesting, because clas- 
sical theoretical work predicts that, for constant mass-to-light ratio, 
r^' " ligh t-profiles have velocity dispersion min ima in the centres of 
galaxies ( lBinnevll98CI : ICiotti & Lanzonil 19971), un less they contain 
central black holes (e.g. iMerritt & Ouinlarj|l998l) . The central a- 
drops evidently do not occur frequently in real early-type galaxies, 
but are, perhaps marginally more common to fast rotators. About 
10% of the SAURON early-types exhibit the central drop (but addi- 
tional 20% have flat central profiles), which is much less compared 
to 46% among the Sa bulg es, also observed with SAURON data 
JFalcon-Barroso et aklboOo) . 



4.3 Distribution of axial ratios 

We now compare average photometric and kinematic axial ratios of 
SAURON galaxies. The axial ratio of a velocity map is related to 
the opening angle of the iso-velocity contours. In other words, the 
pinching of the contours in a spider diagram is related to the axial 



ratio of the best fitting ellipse given by kinemetry. As the kinematic 
axial ratio of slow rotators is an ill defined quantity (set to 1), in the 
rest of this section we focus on the average axial ratios of the fast 
rotators. 

In Fig.|4]we compare values of (qtin) and {qphot} for fast ro- 
tating galaxies. Since the typical seeing for SAURON data ranges 
up to 2'.'5, we exclude the inner 5" of the qtin profiles from our 
derivation of the luminosity weighted average values (AppendixlAt. 
The left hand panel shows a one-to-one correlation between the two 
quantities, although the scatter and uncertainties are large. The right 
hand panel shows more clearly the amount of scatter in these rela- 
tions, as measured by the difference {qtin) - (qphot). The typical 
variation of the measured average Aqtin is 0.1, as shown with verti- 
cal guidelines on the right panel. Outside this region there are about 
dozen galaxies. A few of these have (qkin) larger than (qphot) ', their 
photometric axial ratio is flatter than the kinematic, while the ma- 
jority of outliers have the kinematic axial ratio flatter than the pho- 
tometric. Let us consider in more detail only objects at significant 
distances from the vertical lines (i.e. (qk) — (qp) > 0.15): NGC82I, 
NGC4270, NGC4473, NGC462I, NGC5308 and NGC5838. 

If we look at the qkin profiles of the three galaxies (Fig. IBlt 
with (qkin) > (qphot), we can see that some of their kinematic 
sub-components have axial ratios very similar to the local photo- 
metric axial ratios. Notably, in the case of NGC5308 this is the outer 
component, especially at radii near to the edge of the SAURON 
FoV. The middle range, where the differences are the largest, is also 
the region of the transition between the two kinematic components. 
The mixing of the components changes the measured (qkin), but it 
should be also noted that the (qphot) varies over the whole map, be- 
coming flatter and similar to (qkin) with radius. This is not the case 
in NGC4473. Here the photometric axial ratio remains constant, but 
the kinematic axial ratio changes at larger radii. Again this change 
occurs in the transition region between the two kinematic compo- 
nents. Detailed dynamical modelling of this galaxy shows that it is 
made of two counter-rotating stellar components of unequal mass. 
This object is physically similar to NGC4550, where the main dif- 
ference is in the mass fraction of the two components (Paper X). In 
the case of NGC4270 it is the qphot that steadily changes with ra- 
dius, while qkin has mostly high values, but shows abrupt changes 
between the points. These are related to sometimes rather high val- 
ues of ks/ki, which also changes abruptly between the adjacent 
points, a behaviour originating from the noisy transition region be- 
tween two kinematic components. Since it is hard to disentangle the 
noise from the genuine physical signal in qkin measurement, and 
given the boxy appearance on the large scales, we note that this ob- 
ject could be a true and unusual outlier from the relation between 
qphot and qkin. 

Galaxies with (qkin) < (qphot) are either single component 
(NGC821 and NGC462I, if we exclude the CRC in NGC462I of 
~ 4" in size) or multi component (NGC5838). The HST image of 
NGC5838 has a prominent nuclear dust disk with the axial ratio 
between 0.3 ~ 0.4, constraining its inclination to about 70°. It is 
possible that the velocity map is dominated by the presence of an 
associated stellar disk embedded in the galaxy body. NGC82I and 
NGC4621 were parametrised via the Multi-Gaussian Expansion 
JEmsellem etai]|l994lCapDellarill2002l) in Paper IV. Both models 
required very flat Gaussians to reconstruct the light distribution (in 
both cases, the smallest axial ratios of the Gaussians of the MGE 
models were 0.3); these Gaussians are tracing disks embedded in 
spheroids. If we compare the (qkin) with the flattest MGE Gaus- 
sians, both galaxies move well within the two vertical lines on right 
hand panel of Fig. |4] (blue lines). 
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Figure 4. Left: Relation between luminosity-weighted kinematic and photometric axial ratios, (qkin) and (qphot). respectively, for fast rotating galaxies in 
the SAURON survey measured beyond 5" to avoid seeing effects. The black line is the 1:1 relation. The error bars describe the radial variations of qphot and 
q J. J „ profiles. Right: Difference between luminosity-weighted kinematic and photometric axial ratios plotted against the photometric axial ratio. Dashed vertical 
lines represent a typical variation of qkin profiles, which is the dominant uncertainty factor for comparison with qphot- Solid symbols are fast rotating galaxies, 
while open squares are KDC components of slow rotators. The red lines between symbols link (qkin) with (qphot) measured on inner (NGC4473) and outer 
(NGC5308) kinematic components. The blue lines link (qkin) and axial ratio of the MGE models. 



As a matter of interest, on the same right hand plot of Fig.|4]we 
overplotted axial ratios for the big KDCs in the sample: NGC3414, 
NGC3608, NGC5813 and NGC5831. The shape of their kinemat- 
ics is similar with the distribution of light, except in the case of 
NGC3608. Its KDC has flatter kinematics than the light, but the 
flattest MGE Gaussian is comparable with the kinematic axial ra- 
tio. This suggest that even some of the sub-components of the slow 
rotating galaxies have similar properties like fast rotating galaxies. 

Having in mind that qkin profiles can vary significantly with 
the radius, and the average values could be contaminated by the 
contribution of the transition regions between the components, we 
conclude that there is a near one-to-one correlation between aver- 
age kinematic and photometric axial ratios in fast rotating galaxies, 
with a number of objects having (qkin) < (qphot), and, hence, 
having disk-like components more visible in their kinematics than 
in photometry. 



Shape differences between velocity dispersion and 
surface brightness maps 



4.4 



As shown in Section [Jm the isophotes are not necessary a good 
representation of contours of constant velocity dispersion and the 
deviations are visible in the second cosine term (62) of the har- 
monic decomposition of the velocity dispersion profiles extracted 
along the isophotes. On Fig.|5]we quantify the differences between 
isophotes and iso-cr contours by plotting the normalised luminosity- 
weighted second term ((62/0.0)0-) extracted along the isophotes of 
SAURON galaxies. Slow rotators are shown in red. 

Focusing on the left hand panel, it is clear that galaxies in the 
sample span a large range of (b2/ao)cr, both positive and negative. 
If we exclude NGC4473 and NGC4550, there is a tail of galaxies 



with negative values up to ~ —0.1. Also, it seems that there is a 
trend of high negative {62 /ao ) o- in flat galaxies: as galaxies become 
rounder, (62/00)0- tends to zero; above (qphot) ~ 0.8 galaxies have 
small absolute value of (62/00)0-- Slow rotators are relatively round 
systems (Paper X) with qphot > 0.7 (excluding the special case of 
NGC4550). Their (62/00)0 values are small and mostly positive. 
Typical measurement error of (62/00)0 is ~ 0.03, and hence, the 
slow rotators are consistent with having velocity dispersion maps 
very similar to the the distribution of light, with possibly marginally 
flatter iso-cr contours. 

Before we turn to the right panel on Fig. [5] let us go back to 
Fig. [Hand the example of NGC2549 and NGC4473. Their stellar 
velocity dispersion maps have different shapes from the distribu- 
tions of light, but they have similar absolute values of (62/00)0; 
they fall on the opposite sides of Fig. [5] There is, however, evidence 
that in both cases the deviations from the photometry have similar 
physical origin. NGC2549 shows clear photometric and kinematic 
evidence for a stellar disk viewed at a high inclination (Fig. IB It . 
The velocity dispersion map has a very specific 'bow-tie' shape in 
the central 10", while outside its amplitude is low everywhere. The 
'bow-tie' shape can be explained with a light-dominating dynam- 
ically cold stellar disk which decreases the observed velocity dis- 
persion. Outside of the disk, the bulge dominates and the observed 
velocity dispersion is again higher. It should be noted that the view 
of the galaxy at a high inclination is crucial for this effect to be seen. 

NGC4473 was discussed in the previous section where we 
stated that it actually contains two sub-components. One of them 
is flat and counter-rotates with respect to the main galaxy body. 
This can explain the rise of the velocity dispersion along the ma- 
jor axis, specifically at larger radii where the contribution of the flat 
component becomes similar to the main body. This is also visible in 
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Figure 5. Left: Relation between isophotal axial ratio and the luminosity weighted average normalised second term ((62/00)0-) of the Fourier decomposition 
of the velocity dispersion profiles extracted along the isophotes. Red symbols are slow rotators. An average uncertainty, describing the radial variation of 62/00 
profiles, is shown in the upper left comer. Large negative (62/10)0- values are typical for iso-cr contours that are rounder than the isophotes, while large positive 
(62/^0)0- values are typical for flatter iso-cr contours than the isophotes. Right: same as left, but now absolute values of (62/00)0 are plotted. In both images, 
NGC4550 was shifted to the left for 0.1 for presentation purposes, as shown by the arrows. 



the mean velocity map, where beyond 10" the velocity starts drop- 
ping. The only significant difference between these two examples 
is in the sense of rotation of the flat sub-components. Co-rotating 
disks or flattened components, viewed at favourable angles will 
likely show negative {b2/ao)cr, while counter-rotating components 
will contribute to positive {62/ao)o- In other words, co-rotation in- 
creases, while counter-rotation decreases the flattening of the iso-a 
contours. 

A similar case is NGC4550, the most extreme outlier, which 
also has two counter-rotating disks. NGC4150 and NGC3032 can 
be put in the same group: the OASIS measurements resolve their 
CRC components. It is interesting to note that the strongest cr-drop 
galaxies (NGC2768 and NGC4382) are both found on the posi- 
tive side of the left hand panel of Fig. [5] admittedly with small 
(62/00)0- (< 0.02), indicating that the shapes of the velocity disper- 
sion maps of (7-drop objects are somewhat different from the shapes 
of galaxies with central plateaus in velocity dispersion profiles or 
bow-tie shape of velocity dispersion maps. The shapes of a"-drops 
are more similar to the isophotes than the shape of iso-cr contours in 
the objects with strong near-to edge-on disks such as NGC2549 or 
NGC3377. If cr-drops, bow-ties and central plateaus are all caused 
by disks, than the difference between them might originate from 
the orientation of the galaxy, where bow-ties are seen at edge-on 
and (T-drops at face on orientations. 

Given the uncertainties, it is hard to argue either way about 
other galaxies with slightly positive (62/00)0 values. Other effects, 
such as the actual size of the components or the presence of dust, 
may play a role on that level. This should be noted, but we, however, 
continue by suggesting that the main contributor to the shape differ- 
ence between the zeroth and the second moment of the LOSVD in 
early-type galaxies are the embedded disks or flattened fast rotating 
components seen at different inclination. In this respect, on the right 



hand panel of Fig.|5]we plot the absolute value of (62/t[o)o- Here a 
trend is clear: from round slow rotators with generally small abso- 
lute values to flat fast rotating galaxies with increasing |(62/ao)o|- 
This sequence is the one of increasing contribution of the ordered 
rotation towards the total energy budget, but it is, unfortunately, de- 
pendent on the viewing angles (NGC524 with the face-on disk in 
the top left corner of the diagram) and on the luminosity-weighted 
contribution of the flat component to the total kinematics (CRC in 
NGC4621 is not detected). 



4.5 Higher order moments of the LOSVD 

The final part of our analysis of the kinematic moments of the 
LOSVD is devoted to maps of hs and /14, the Gauss-Hermite mo- 
ments which measure deviations of the LOSVD from a pure Gaus- 
sian. They were introduced because LOSVD profiles are rather non- 
Gaussian, specifically the contribution of hs, the skewness, along 
the major-axis of galaxies was found to have large amplitudes (e.g. 
van der Marel et al. 1994) . In a diagram of local /13 — V/a relation. 
Bender et alj ( 119941) plotted the major axis points of their objects 
which could be separated in two distinct groups made of galaxies 
with disky isophotes and non-disky isophotes. Both groups showed 
anti-correlation between h^ and V/a, where the points with small 
V/a values in the centre of the diagram had a steeper slope. It was 
recognised that the disky early-type galaxies cover a large range in 
V/a and are responsible for tails in upper left and lower right comer 
of the diagram, while non-disky early-types with small V/a values 
showed a somewhat smaller range in /13 values. 

On the top panel of Fig. |6] we plot the local /13 — V/a rela- 
tion for all data points of 48 S AURON galaxies. Red points belong 
to Voronoi bins of the slow rotators. This plot is somewhat differ- 
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Figure 6. Top: Local /13 — y/cr relation for all data points of the 48 
SAURON galaxies with h^ eiTor less than 0.2. Middle: Profiles of ^3 - V/a 
relation averaged on ellipses. On all maps (V , cr and ^3) profiles were ex- 
tracted along the best fitting ellipses to the velocity maps. Typical en'or bars 
are of the order of 0.01 or smaller, due to averaging along the best fitting 
ellipse. Bottom: Relation between the luminosity averaged /13 and {V/cr)e 
(V/cr within 1 Re from Paper X). The error bars des cribe the radial varia - 
tion of /13 profiles. The solid line is the relation from lBender et alj 119941) . 
On all panels, slow rotators are sho wn in red. 
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ent from previous findings. Paper X showed that slow rotators have 
small global V/a and this is also reflected in bin by bin values. The 
black points represent the Voronoi bins of fast rotators, and their 
distribution is different from what was found before for disky ellip- 
ticals. The shape of the distribution of black points can be described 
as a superposition of two components: one which is anti-correlated 
with V/a and makes distinct tails in upper left and lower right cor- 
ner of the diagram (large positive and negative V/a values), and 
the other that is both consistent with hs ~ 0, and shows positive 
correlation at intermediate V/a. 

Existence of /13 ~ distribut i on of points could be explained 
by the fact that, while lBender et alj ( 119941) plot only the points along 
the major axis, we plot them all, and the horizontal distribution and 
correlating tails are the consequence of the minor axis contamina- 
tion. Using kinemetry we extracted profiles of the dominant har- 
monic terms from V, a, and ha maps along the same best-fitting 
ellipses to the velocity mapo We plot the ha — V/a profiles on 
the middle panel of Fig. [6] These profiles are basically major-axis 
representation of two-dimensional maps and are more comparable 
with previous major axis data. The slow rotators (in red) have the 
smallest V/a values and generally small amplitudes of ha . The two 
trends for fast rotators from the left-hand panel are still visible. Fast 
rotators cover a range of V/a values, but some have large and some 
small ha amplitudes for a large V/a values. In some specific cases 
there is a suggestion that ha changes the sign becoming positive and 
correlating with V/a (e.g. NGC5308) at larger radii. 

A confirmation of this can also be seen on the bottom panel of 
Fig. [6] Here we plot the relation between luminosi ty averaged ha 
and {V/a)e obtained within I Re, as advocated bv Isiimevl ( 120051 
values taken from Paper X). The shown error bars are not statitcial 
uncertainties but describe the radi al variation of ha p rofiles. This 
plot can be compared to Fig. 14b of lBender et alJ ( 119941 ). which also 
shows galaxies with {/13) ~ for intermediate V/a^, while their 
empirical fitting relation (the solid line) describes the general trend 
in our data. 

The last kinematic moment to be analysed is given by the /14 
Gauss-Hermite coefficient. This moment describes the symmetric 
departures of the LOSVD from a Gaussian. It should be, however, 
kept in mind that it is very difficult to measure /14 robustly, since 
it strongly depends on the effect of template mismatch. Also an 
inaccurate removal of the continuum will cause spurious h^ values. 

On the left hand panel of Fig. |7]we plot the local /i4 — V/a 
relation. Slow rotators with small V/a values, dominate the central 
region. They cover a range of positive h^ values and somewhat ex- 
tend below zero. Fast rotators fill in a cloud around slow rotators, 
equally filling negative and positive V/a part of the plot. There is a 
suggestion for a larger spread in V/a for positive /14 values. Lumi- 
nosity weighted average (hi) values extracted with kinemetry along 
the isophotes are presented as a function of luminosity-weighted 
V/a measured within 1 Re (Paper X) in right panel of Fig. |7] As 
{V/a)e increases, there is a marginal trend of an increased spread 
in the observed /14 values, and galaxies wi th negative average val- 
ues start to appear. This result is similar to lSender et alj jl994h . in 
the sense that negative (h^) appear for larger {V/a)e, but we do not 
find as negative values of (/14). The only slow rotator with negative 



^ Note that here we re-run kinemetry on a maps along the best fitting el- 
lipses for velocity maps, except in the case of the slow rotators when we 
used circles as explained in Section [T2l These profiles are somewhat differ- 
ent from the profiles presented in Section l4!4l which were obtained along the 
isophotes. 
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Figure 7. Right: Local correlation between /14 and V/cr for all data points of the 48 SAURON galaxies with /14 error less than 0.2. Right: Correlation between 
luminosity weighted values average (h^) and {y/c!)e (Paper X). Red symbols are slow rotators in both panels. The en'or bars describe the radial variation of 
/i4 profiles. 



(hi), however, is the usual outlier, NGC4550, the special case of 
two counter-rotating disks. 



5 DISCUSSION 

5.1 Velocity maps of slow and fast rotators 

Radial profiles of kinemetric coefficients show that early-type 
galaxies are (i) multicomponent systems and (ii) in the majority 
of the cases contain a kinematic equivalent of a disk-like compo- 
nent. These statements are based on the empirical verification that 
the assumption of kinemetry holds. The assumption is that the az- 
imuthal profiles, extracted from a velocity map along the best fit- 
ting ellipses, can be described with a simple cosine variation. Prac- 
tically, this means that higher-order harmonic terms are negligible 
and at our resolution we find that the deviation from the pure co- 
sine law is less than 2%, for about 80% of cases, at least on a part 
of the map. The multiple components are visible in the abrupt and 
localised changes of kinemetric coefficients. The central regions of- 
ten harbour separate kinematic components which co- or counter- 
rotate with respect to the outer body. In some cases there are more 
than two components (e.g. NGC4382) or components are similar in 
size and co-spatial (e.g. NGC4473). Although the axial ratio pro- 
file in edge-on systems can change with the seeing effects (see Ap- 
pendix [Ajl, and in some cases the seeing can alter the map consid- 
erably, the kinematic sub-components are usually robust features. 

There are, however, early-type galaxies for which the devia- 
tions from the cosine law exceed 10% across a significant radial 
range. Such objects are all classified as slow rotators in Paper IX. 
The breakdown of the kinemetry assumptions is another evidence 
that these objects are intrinsically different. Certainly, we might not 
be able to apply kinemetry a priori in its odd version on veloc- 
ity maps that do not show odd parity (e.g. NGC4486), but in the 



case of slow rotators that show a detectable level of rotation (e.g. 
NGC5982), the kinemetric analysis clearly shows differences from 
the maps of fast rotators (e.g. NGC4387). 

There are three types of intrinsic structures that will show 
small mean velocities: (i) face-on disks, (ii) two counter-rotating 
equal in mass and co-spatial components and (Hi) triaxial structures. 
(i) Face-on Disks. Low inclination thin disks will still have minimal 
deviations from the cosine law assumed in kinemetry. Although the 
amplitude of the rotation in such cases is small, the disk-like ro- 
tation (negligible ks/ki coefficients) is expected at all projections 
except in the actual case of i = 0° with no rotation. Such an exam- 
ple in the sample is NGC524 (i ~ 19°). Based only on the velocity 
map, one could naively interpret that NGC4486 is a thin disk at 
i = 0°, given that it has qtin ~ 1 and shows no rotation. How- 
ever, other moments of the LOSVD strongly rule out this geometry, 
while, in general, other slow-rotators are not round enough to even 
be considered as such extreme cases. 

(ii) Counter-rotating components. This case can be outlined with 
NGC4550 and NGC4473. These galaxies are examples of axisym- 
metric objects with two counter-rotating and co-spatial compo- 
nents. The first one is classified as a slow rotator, while the other 
one is a fast rotator. In both cases, however, it is the mass fraction of 
the components that really decides what is observed. In the case of 
NGC4550 the masses are nearly equal and the luminosity weighted 
mean velocity is almost zero in the central region. NGC4550 is a 
product of a very specific formation process, but a clear example 
of how a superposition of two fast rotating components can imi- 
tate a slow rotator. In NGC4473, on the other hand, one compo- 
nent is more massive and dominates the light in the central 10", 
where the rotation is clearly disk like. Outside this region, where 
the counter-rotating component starts to significantly contribute to 
the total light, the shape of the velocity map changes, the amplitude 
of the rotation drops and non-zero kz,/k\ coefficients are necessary 
to describe stellar motions. 
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Figure 8. Left: Difference between luminosity weighted kinematic and photometric axial ratio, {qkin) and {qphot)^ for 5 Jeans models at 25°, 35°, 50°, 
65°, 80° and 90° inclinations plotted against the luminosity weighted photometric axial ratio. The increase in inclination decreases {qphot) and moves points 
downwards. The vertical dashed lines are guidelines of the variation of qkin profiles from the SAURON data as in Fig.|4] Vertical solid line shows the location 
of zero difference between (qphot) and {qtin)- Right: Relation between photometric axial ratio and absolute (62/00)0- for the same Jeans models as on the 
left. With increasing inclination the points move towards more negative (b2/ao)„- values. The hatched area represents the location of the data from Fig.|5] 



(Hi) Triaxial objects. Observationally, these objects are marked 
by kinematic twists and kinematic misalignments, which are not 
present in axisymmetric galaxies. A strong restriction to the shape 
of the velocity maps of axisymmetric galaxies is that they should 
not have a radial variation of PAkin or a misalignment between 
kinematic and photometric position angles. This, however, is not 
the case for triaxial galaxies, where the change of position angle 
can be influenced by a true change of the angular momentum vec- 
tor, by the orientation of the viewin g angles or by the relative dom- 
inance of different or bital families dFranx et al Jl 1 99 ll : Istatleill 1 99 it 
Ivan den Bosch et al.t2008.) . Paper X found that, globally, kinematic 
misalignments are present only in slow rotating galaxies. Locally, 
however, we find both kinematic misalignments and twists in some 
fast rotating objects as well, but they are mostly confined to central 
components or clearly related to bars (e.g. NGC1023) or galaxies 
with shells (e.g. NGC474). 

We do not find strong kinematic twists typical of the extreme 
cases of maximum entropy models of tri axial galax ies projected 
at var ious viewing angles (( IStatlenll99ll) ; see also lAmold et alj 
( Il994h ) for velocity maps of Stackel triaxial models). Given that 
our sample is not representative of the luminosity function of lo- 
cal early-type galaxies in the sense that it contains too many mas- 
sive galaxies, which are also more likely to show extreme features 
on the velocity maps, it is remarkable that we find that only a few 
velocity maps are similar to those predicted. Still, the observed ve- 
locity maps are divers (e.g. the difference between the maps of slow 
and fast rotators) and their complexity reflects the difference in their 
internal structure. 

The profiles of the relative change of the kinematic position 
angles from the first panel of Fig. [3] suggest that fast and slow ro- 
tators have genuinely different intrinsic shapes, fast rotators being 
mostly axisymmetric and slow rotators weakly triaxial. This is also 



reflected in k^/ki ratio. We suggest that the high values of this ratio 
in slow rotators, which is in practice caused by the noise in non- 
rotating velocity maps, has its origin in the internal orbital make up 
of these galaxies. Weakly triaxial slow rotators contain box orbits, 
and competing contributions of different tube orbit families, as op- 
posed to more axisymmetric fast rotators wit h short-axis tubes as 
the only major orbit family jde Zeeuwl 1 198511 % This suggestion is 
also supported by the analysis of the orbital structure of collision- 
less merger remnants (Jesseit et al. 2005) as well as by the kine- 
metric an alysis of velocity m aps of simulated binary disc merger 
remnants jjesseit et alj2007l) . 



5.2 Evidence for disks in fast rotators 

Kinematic sub-components with azimuthal profiles that can be fit- 
ted with a cosine law are described as having a Disk-like rotation 
(DR). This does not mean that they are actual disks. It just sug- 
gests that velocity profiles of early-type galaxies extracted along 
the best-fitting ellipse resemble the velocity maps of thin disks in 
circular motion. The rate of occurrence of DRs is, however, strik- 
ing. There is no reason why this should be the case in early-type 
galaxies, which in principle as a class can have a triaxial symme- 
try and complex motion in different planes. As suggested above, 
the link between the kinemetry assumption and the structure of fast 
rotating galaxies has its origin in their internal structure. 

The results of the dynamical models in Paper X reveal 
that fast rotators show evidence for a kinematically distinct flat- 
tened spheroidal component, suggesting that fast rotators are 

* We should keep in mind that two short-axis tube families with oppo- 
site angular momentum in certain cases can produce axisymmetric objects 
which appear as slow rotators, as mentioned above in the case of NGC4550. 
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nearly oblate and contain flattened components. In addition to 
these dynamically cold components, the stellar populations of 
fast rotators show evid ence for different chemical components. 
iKuntschner et alj ( 120061) find that all morphologically flat fast ro- 
tators have Mgb line-strength distribution flatter than the isophotes, 
and associate it with the rotationally supported substructure, which 
features a higher metallicity and/or an increased Mg/Fe ratio as 
compared to the galaxy as a whole. 

These are some of the dynamical and chemical evidences for 
disk-like components in fast rotators. What is the kinemetric evi- 
dence? As mentioned above, we find that velocity maps of fast ro- 
tators are mostly described by a simple cosine law, as are velocity 
maps of thin disks. We also find an almost one-to-one correspon- 
dence between the projected shape of the stellar distribution and the 
shape of the observed kinematic structure in fast rotating galaxies. 
The connection between the shape and the kinematics is supported 
by an assumption that rotation influences the shape of the object by 
flattening it and, for isotropic models, the rotation speed responsi- 
ble for flattening of t he shape is related to the s hape of the stellar 
distribution as ~ ^ JBinnev & Tremainell2008l . Section 4.8.2). In 
order to investigate further the (qphot) — (qfcm) correlation we con- 
structed two-integral analytic models of early-type galaxies. 

The isotropic models we used were previously presented in 
Appendix B of Paper X, to which we refer the reader for more de- 
tails. The main point of these Jeans models is that we used as tem- 
plates 6 galaxies, which represent some of the typical types from 
the SAURON sample. Their light distribution was parameterised 
in Paper IV by the MGE method, and was used as the basis for 
the intrinsic density distributions. Observables of each Jeans model 
were projected at 6 different inclinations: 90°, 80°, 65°, 50°, 35° 
and 25°. These models are not meant to reproduce the observed 
kinematics in detail, but they are self-consistent, and under the as- 
sumption of axisymmetry and isotropy, they predict velocity maps 
and offer an opportunity to study the relation between the shape and 
kinematics. 

On the left hand panel of Fig. [8] we show the difference 
between the luminosity-weighted average values of (qfcin) and 
(qphot), measured by kinemetry on the model images and veloc- 
ity maps, in the same way as for the SAURON data in Fig. |4] (also 
excluding the inner 5"). Different colours represent Jeans models 
based on different template galaxies. Each symbol corresponds to 
a model at different inclination, where the points move from top 
to bottom with increasing inclination (from 25 — 90°). It is clear 
that for small inclination {qfci„) « (qphot), but as the models are 
viewed closer to edge-on there is a trend of increasing differences 
between the axial ratios and in some cases a trend of larger variation 
along the profiles represented by larger error bars for progressively 
more inclined models. Specifically, in all but one marginal case 
(qfcin) > (qphot): in these models velocity maps are 'rounder' than 
images. The velocity map of the Jeans model of NGC4621, whose 
MGE parametrisation has the flattest Gaussian, has a flat component 
along the major axis, which becomes more prominent with increas- 
ing inclination, contributing to the radial variation and increasing 
qfcm with respect to qphot- The kinematics of this component is 
the most 'disk-like' in our Jeans models with tightly pinched iso- 
velocity contours. For comparison, the model for NGC3377 also 
has a similar thin MGE component. This component contributes to 
the disk-like kinematic component confined to the central region, 
but since it is not as prominent in the total light, such as the one in 
the NGC4621 MGE model, only traces of the disk-like component 
can be seen in both photometry and kinematics. 

The contrast between observed galaxies and isotropic models 



is significant. The isotropic models predict (q^in) > (qphoi) for 
i > 30°, but the models also show that prominent disk-like pho- 
tometric features will generate pinched iso-velocity contours and 
decrease (qfein) pushing the galaxies towards the observed trend of, 
on average, {quin) ~ {({phot) or even {qfci„) < {({phot)- It seems 
reasonable to assume that while isotropic models can explain cer- 
tain features of fast rotators, they cannot describe them as a class 
of galaxies. It is the embedded flattened components, often visible 
only in the kinematics that are responsible for the observed differ- 
ences between the models and the data. 

This can be also seen by comparing the difference in shape 
between the isophotes and iso-cr contours (Fig. |5}. We repeated 
the same exercise with velocity dispersion maps of our isotropic 
Jeans models. Result are shown on the right panel of Fig. [8] The 
trend shown here of a larger absolute {b2/ao)<y values with lower 
{({phot) is very similar to the observed trend (hatched region). The 
test can also explain the shape of the observed trend: large abso- 
lute values of {&2/ao)<T can be observed when the object contains a 
significant flattened component and it is observed at larger inclina- 
tions. The most affected are again NGC4621 and NGC3377 mod- 
els. As before, the isotropic models are able to explain part of the 
observed data, but not the details of the distribution. Specifically, 
the shape of the iso-o" contours of slow rotators {{h2/ao)<j > 0) is 
not reproduced well by the isotropic models. Similarly, the spread 
in (62/10)0- of fast rotators is also not well reproduced. Clearly, 
our Jeans models are much simpler than the real galaxies, lacking 
by construction multiple kinematic and especially counter-rotating 
components. Comparing the kinemetric analysis of the Jeans mod- 
els and the observed objects, we find that fast rotating galaxies 
are more complex than isotropic rotators, presumably containing 
also flattened kinematically distinct components, which can co- or 
counter-rotate on top of the non-rotating or isotropically rotating 
spheroid. 

The evidence for disks in fast rotators are also present in the ra- 
tio of /i3 and V/o. hz measures the asymmetric deviations from the 
Gaussian LOSVD and the anti-correlation of /13 with V/a is taken 
to show presence of disks. Our results confirm previous findings 
that early-type galaxies on the whole have asymmetric LOSVDs, 
but this applies only to fast rotators. We also find that many fast 
rotators show constant and close to zero hz profiles, or, in a few 
cases, show a change from negative to positive values with in- 
creasing r adius, similar to what is se en in peanut bulges and bars 
JChung & Bureau 2004 ; Bureau & Athanassoula 2005) 

The embedded flattened components in fast rotators are also 
evident when /13 — V/o diagram is compared with the re- 
sults of merger simula t ions (iBalcells: 1991: Be ndo & Bamesll2000l ; 
iNaab & Burker3l200ll : iGonzalez-Garcia et al.l 12006). Specifically, 
the updated hz — V/a diagra m in Fig. |6] can be compared with 
Fig. 16 from lNaab et alj ( 120061 ) ■ who discuss the influence of dis- 
sipational mergers in which embedded disks are formed in merger 
remnants. Our figure compares rather well with a combination of 
1:1 dry and 1:3 wet mergers. This comparison suggest that it is 
not possible to explain the LOSVD of early-type galaxies with one 
merging track only, but that slow rotators predominantly originate 
in major colissionless mergers, while fast rotators are remnants of 
dissipational mergers. 

The comparison of our bin-by-bin /14 — V/a dia gram (right- 
hand panel in Fig.|7} and lower four panels of Fig. 16 in lNaab et al.l 
( 12006.) is equally impressive, although their simulations predict 
somewhat too negative values of /14. Again, the observations are 
largely consistent with the scenario where slow rotators originate 



© 2007 RAS, MNRAS 000,[TH2| 



The SAURON project Xn 17 



from dry 1 : 1 mergers, while fast rotators from a combination of dry 
and wet 1:3 mergers. 

Combining all the evidence presented in the previous sec- 
tion, we suggest that fast rotators are dominated by disks (e.g. 
NGC3156, NGC2685). When their light is dominated by the bulge, 
their kinematics still show strong disk components (e.g. NGC821, 
NGC4660). In either case, fast rotators contain flattened fast ro- 
tating components and this dynamical property differentiates them 
from slow rotators. We suggest that with increasing specific angular 
momentum, Ajj, the relative mass of the embedded disks also in- 
creases and contributes more significantly to the total mass. Among 
the disk-like components in fast rotators there is a range of flatten- 
ings reflecting a diversity in possible formation paths which create, 
preserve and/or thicken disks within spheroids, such as passive fad- 
ing of spirals and multiple minor dissipational mergers. The change 
within the internal structure is observationally reflected in the tran- 
sition between slow and fast rotators, which offers possible anchor 
points for theoretical models of galaxy evolution. 



6 CONCLUSIONS 

Using kinemetry we analysed two-dimensional maps of 48 early- 
type galaxies observed with SAURON and velocity maps of a 
sub-sample observed with OASIS. The analysed maps are: re- 
constructed image, mean velocity, velocity dispersion, h^ and /14 
Gauss-Hermite moments. The reconstructed images and the maps 
of the mean velocity were analysed along the best fit ellipses. Maps 
of the higher moments of the LOSVD were analysed either along 
the isophotes (velocity dispersion, hi) or best-fitting ellipses from 
velocity maps (/13). We presented the profiles of kinemetric coeffi- 
cients for velocity maps, being the dominant kinematic moment of 
the LOSVD and having the highest signal-to-nose ratio. 

Kinemetry and its kinemetric coefficients allow us (i) to dif- 
ferentiate between slow and fast rotators from velocity maps only 
and (ii) to indicate that fast rotating galaxies contain disks with a 
larger range in the mass fractions to the main body. The results of 
this work can be summarised as follows: 

• Following the kinemetric analysis of the velocity maps it is 
possible to distinguish between galaxies with Single and Multiple 
Components. Components can be described as having a Disk-like 
rotation. Low-level velocity. Kinematic misalignment. Kinematic 
twist or being Kinematically Decoupled Component. The sorting of 
galaxies in these groups is based on kinemetric coefficients: PAfci„, 
qfc,„, fci andfcs/fci. 

• The majority of early-type galaxies are MC systems (69%) 

• The total fraction of galaxies with a DR component (including 
afl that have fcs/fci < 0.02) is 81%. In terms of SO/E classifica- 
tion, 92% of SOs and 74% of Es have components with disk-like 
kinematics. 

• KDCs are found in 29% of galaxies. These KDC are of differ- 
ent sizes, some are not resolved at SAURON resolution and appear 
as CLVs, but are clearly detectable in the OASIS observations at 
higher spatial resolution. 

• Early-type galaxies with constant PAfei„(APA < 10° outside 
any central component) are fast rotators. All fast-rotators have a DR 
component. 

• Most fast rotators have {(\kin) ^ (qphot) • Their kinematics of- 
ten show a structure flatter than the distribution of light. This means 
that images alone are not sufficient to recognise all fast rotators (e.g. 
NGC524 or N3379 would be missed) 



• In face-on galaxies, isophotes are coincident with iso-a con- 
tours. In edge-on galaxies, however, there are differences which can 
be detected with kinemetry. Specifically, the edge-on fast rotators 
contain dynamically cold components which changes the shape of 
iso-a contours. 

• Slow rotating galaxies have low V/a and a range in h^ am- 
plitudes. Fast rotating galaxies have large spread in /13 which anti- 
correlates with V/a. There is, however, a significant number of fast 
rotating galaxies which have h-j, radial profiles that change from 
strongly negative to near to zero /13 values. Similar behaviour was 
found in bar galaxies or in remnants of colissionless mergers. 

• Allowing for large uncertainties in the determination of hi, 
our data show a trend where slow and fast rotators have similar 
hi values, with a weak tendency for an increased spread of hi in 
fast rotators, also dependant on V/a. These trends can be explained 
through a combination of dry and wet mergers of both equal and 
unequal in mass progenitors. 

• Dissimilarities between slow and fast rotators originate in their 
different internal structures. Slow rotators are mildly triaxial objects 
supporting a variety of orbital families. Fast rotators are axisym- 
metric spheroids with embedded flattened components of different 
mass fraction ranging from completely disk dominated systems to 
small central disks visible only in kinematics. 
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APPENDIX A: INFLUENCE OF SEEING ON THE SHAPE 
OF THE VELOCITY MAPS 

The effects of seeing on galaxy images in the regime of non- 
adaptive optic as sisted observations are now well understood 
JSaglia et al.ll 19971) . They are specifically important in the nuclear 
regions where the instrumental effects and the atmospheric see- 
ing redistribute photons sufficiently enough to change the intrinsic 
ellipticity or position angle an d mask structures like concentrated 
stella r nuclei or do uble nuclei JSchweizegl 19791 . 1 198 u \ iDiorgovskil 
1 19831 : iLauen 119851) . Specifically, luminosity profiles are affected 
more by the seeing when the ellipticity is low, while the ellipticity 
profil es are affected more when the ellipticity is high jPeletier et alj 
1 1990) . In general, effects of seeing can be detected to quite large 
radii. 

It is natural to expect the redistribution of photons due to the 
seeing will also influence the two-dimensional shape (maps) of the 
higher moments of the LOS VD. Specifically, for the first moment, 
the mean velocity, the obvious effect of larger seeing is a less steep 
rise of the velocity in the centre. Additionally, one can expect an 
increased axial ratio (rounder map), or in other words, a more open 
spider diagram in the central regions. Except for these qualitative 
expectations there are no quantitative estimates of the influence of 
seeing on the shape of velocity maps. Observations of the SAURON 
sample were taken under a variety of atmospheric conditions, so 
before analysing the velocity maps, we tested the effects of different 
PSFs. 

We construc ted a number of model velocity maps using a 
iHemguislI J 19901) potential which reasonably well approximates 
the density of real early- type galaxies. Each model map was con- 
structed from the Hernquist circular velocity profile, projected at a 
certain inclination to imitate different axial ratios seen in the sam- 
ple (qfcin), and weighted by the given Hernquist surface brightness 
profiles. All models have the same scale length ro = 10". The in- 
clinations ranged from 20 to 75 degrees, where and 90 degrees are 
face on and edge on viewing angles, respectively. This corresponds 
to velocity maps with axial ratios quin in the range of 0.94 to 0.26. 
Although in the SAURON sample there are near to edge on galax- 
ies, their velocity maps do not resemble the limiting case of edge 
on disks (qfci,i=0), being considerably less flat. The above range in 
axial ratios were used to match the values observed in the SAURON 
sample. Maps were then convolved with a Gaussian kernel of dif- 
ferent full-wi dth-half-maximum (FWHM), accounting also for the 
square pixels JOianet al.l 19951) in such a way as to make the model 



maps similar to observations, especially in the central regions which 
remain unbinned. The range of FWHM was also matched to the 
measured seeing range of SAURON observations (Paper III). Fig- 
ure lAll shows an example of model maps for axial ratio of 0.5 and 
seeing FWHMs of 1", 2" and 3". 

We ran kinemetry on the convolved model velocity maps and 
extracted the same parameters as for the observed velocity maps: 
PAfci„, quin, k\ and k^/k\ harmonics. The results can be seen in 
Fig. lA2l and they suggest that: 

• determination of the position angle is not influenced by the 
seeing; 

• the influence of the seeing on the determined axial ratio is 
larger on maps with intrinsically smaller axial ratio. Also, the bigger 
the seeing is, the rounder (larger axial ratio) is the observed map. 

• the influence of the seeing on the rotation curve (fci term) is 
larger on maps with intrinsically smaller axial ratios, where the big- 
ger is the seeing, the lower is the maximum velocity. The position 
of maximum velocity is also pushed towards larger radii. 

• the influence of the seeing on the higher terms (fcs/fci) is negli- 
gible in most cases, except for intrinsically very flat maps, but even 
there the signal is of the order of 1% or less and is not detectable 
with the current accuracy. 

As expected, the seeing influences the axial ratio of velocity 
maps (or, in other words, the opening of the iso- velocity contours), 
and the steepness of the rotation curve, both in the maximum values 
and the positions of these maxima. 

The largest effect is for maps with small axial ratio (large open- 
ing angles), which can be understood as the consequence of the con- 
volution, similar as in the case of photometry and high ellipticity. In 
terms of changing of the intrinsic values, it is the axial ratio that is 
the most affected by seeing. At an intrinsic flattening of 0.77, a bad 
seeing of 2'.' 5, will increase the flattening by 0.1. This effect rises 
to surprising 0.5 difference in measured flattening for an intrinsic 
flattening of 0.26. This effect is milder for the maximum rotation 
velocity, where only for velocity maps with the smallest axial ratios 
and strongest seeings the difference between the intrinsic and see- 
ing convolved values exceeds 10 km/s. Another effect of seeing on 
the velocity maps is the change of position of the maximum axial 
ratio and the maximum rotational velocity with increasing seeing. 
In both cases the trend is clear: larger seeing FWHM causes dis- 
placement of the maximum {qkin or fci) towards larger radii. 

In summary, seeing can strongly influence the appearance of 
the velocity maps, both in the shape, value and extent of the fea- 
tures. This is important for interpreting the structures on the maps 
and, especially, when looking for kinematic subcomponents. 



APPENDIX B: KINEMETRIC PROFILES 

In this section we present the results of kinemetric analysis of 
SAURON velocity maps. 
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Figure Al. An example of model velocity maps for a given inclination and convolved with different seeing kernels. From left to right: A model velocity map 
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Figure A2. Kinemetric analysis of seeing convolved model velocity maps. From left to right: Kinemetry coefficients from maps with axial ratios of 0.94, 
0.77, 0.5 and 0.26 degrees. From top to bottom: Kinemetric coefficients: position angle, axial ratio, fci and k^/ki harmonics. Colours correspond to values 
extracted from maps convolved to different seeing. Black: unconvolved model map. Blue: model map convolved with FWHIVI of 1". Light blue: model map 
convolved with FWHIVI of 1'.'5. Green: model map convolved with FWHIVI of 2". Yellow: model map convolved with FWHM of 2''5. Orange: model map 
convolved with FWHM of 3". Vertical lines show the radial extent of the seeing at FWHM. Note that the incoirect determination of axial ratio at the innermost 
radius is an artifact of poor sampling. 
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